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INTRODUCTION 

There is great need for atmospheric trace constituent measurements with 
higher vertical resolution than attainable with passive radiometers. Infrared 
(IR) DIAL, which depends on Mie scattering from aerosols, has special advantages 
for tropospheric and lower stratospheric applications and has great potential 
importance for measurements from Shuttle (ref. 1) and aircraft. Differential- 
absorption LIDAR data reduction involves comparing large amplitude signals which 
have small differences. The accuracy of the trace constituent concentration 
inferred from DIAL measurements depends strongly on the errors in determining 
the amplitude of the signals. Thus, the commonly used SNR expression (signal 
divided by noise in the absence of signal) is not adequate to describe DIAL 
measurement accuracy and must be replaced by an expression which includes the 
random coherent (speckle) noise within the signal (refs. 2, 3, and 4). A com- 
prehensive DIAL computer algorithm (ref. 5) is modified to include heterodyne 
detection and speckle noise. Results of a parametric study are presented and 
comparisons with direct detection are discussed. Examples are given for mon- 
itoring vertical distributions of , H 2 0, and NH using a ground-, aircraft-, 
or Shuttle-based pulsed tunable CC> 2 laser DIAL system. 

ANALYSIS OF DIAL SENSITIVITY WITH HETERODYNE DETECTION 

The expectation value P of the number of measured photons from one coher- 
ence area of a scattering cell at range R of length Ar is 

P = nrE ArA fi /R -2(C + op) dr (1) 

hVR 2 

where f| = detector efficiency, T = optical efficiency, E = laser energy, 

(3 = 180° backscatter coefficient per length per steradian, Ar = cell length = 

T = integration time, V = frequency, A = transmitter area (= receiver area for 
heterodyne system with single detector) , B, = extinction coefficient (total 
minus that of measured gas), 0 = absorption coefficient of measured gas, and 
p = density of gas being measured. 

In differential absorption, measurements are made of two frequencies 
selected to maximize signal and differential absorption of the species being 
neasured while minimizing interference effects. The double ratio of signals at 
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adjacent scattering cells at two frequencies yields information about absorp- 
tion in the region between the scattering cells. 
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where P^j = expectation value of measured return from cell j at frequency i. 

= backscatter coefficient from cell j at frequency i. Act = absorption 
coefficient difference between frequencies 1 and 2, AC = extinction coeffi- 
cient difference between frequencies 1 and 2 (not including gas being measured) 
i = 1 or 2 frequencies on or off absorption line, respectively, and j = 1 or 2 
for distances R^ and R 2 , respectively. For scattering cells of equal length 
the resolution length R 2 - Rj equals the scattering cell length Ar. The term 
012 021 

In o' ' ~ o ' — ‘ is a correction term due to changes in backscattering with fre- 
Pll P22 R 2 

quency across scattering cells. The term 2/ (AC) dr is a correction due to 

R 1 

interferent species. Variation in backscatter and interferent absorption at 
the two frequencies will result in biases in the inferred concentrations. 

These biases can be reduced by careful frequency selection and partially cor- 
rected by using a priori information and auxiliary measurements. For an opti- 
cally thick species, a series of "on" frequencies is required to maximize sensi 
tivity at various altitudes. 


The random error in p, which is calculated assuming that signal plus 
background are measured during each pulse and that background is measured 
between pulses and subtracted, will depend on the random error in the measure- 
ment of Pf j . The uncertainty in the inferred concentration (6p) 2 is given by 

(3) 


where N is the number of pulse pairs per measurement and SNR is the single 
pulse signal-to-noise ratio. For heterodyne detection, the major errors in 
P^j are due to quantum noise in the local oscillator and fluctuation noise in 
the return signal. 

The heterodyne signal-to-noise ratio for a single coherence volume 

R i j 

is limited to 1.0 due to the speckle noise in the return signal. The 

Pj_j + BT 

number of coherence lengths per scattering cell is BT. The number of coher- 
ence areas viewed by the detector is M, where M is the ratio of receiver to 
transmitter area. For heterodyne detection an individual detector is required 
for each coherence area. The single pulse signal-to-noise ratio for determin- 
ing the signal from a scattering cell is then 

P. . 

SNR = — — - 3 - — ■ /BTM (4) 
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where BTM is the number of statistically independent samples from a scat- 
tering cell for a single pulse. The post-detection bandwidth B is con- 
strained by matching with pulse duration (B < 1/Tp) and by the width of 
atmospheric spectral lines. The integration time T is constrained by verti- 
cal resolution requirements. Neglecting the system error terms, the random 
error in density for heterodyne detection with bandwidth B, integration time 
T, N pulse pairs, and M coherence areas with one detector per coherence 
area is 
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For N pulse pairs and M detectors , the SNR — - — 1 ~ ) - - - - /BTNM is maximized 
^ P i j + BT 

for a given total laser energy per measurement (“ NMPij) when the laser energy 

per pulse per detector is selected so that P^j ~ BT. For that condition an 

approximate solution to equation (5) is 




r yj -L i.'ii'i 


(6) 


This equation illustrates that 6p is independent of p and proportional to 
1/Ao within the constraint that the two-way integrated absorption through the 
atmosphere does not limit Pj_j to less than BT. Equation (6) also indicates 
the strong effect of range resolution on measurement error. The number of 
pulse pairs required to maintain a constant measurement error is proportional 
to the inverse cube of the range resolution. 


COMMENTS ON DIRECT DETECTION 


Direct detection can have advantages over single-detector heterodyne 
detection when signal levels are high since direct detection allows averaging 
over multiple coherence areas with a single detector. The major disadvantage 
of direct detection is background noise which is limited by optical filters 
(n, 10 11 Hz) compared with the electronic filters (10® to 10 9 Hz) for heterodyne 
detection. Reducing the background noise by reducing the field of view will 
result in an increase in speckle noise. 


The direct detection signal-to-noise ratio can be written as a function 
of M, the number of coherence areas, and the total direct detected signal 
p£j = MP-[j. The SNR for direct detection (assuming zero detector noise, 
ref. 3) is 


(snr) dd 
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Equation (7) includes a speckle term Pij /M, a Poisson term P^jBT, and a 
background term (GDMT)BT. In equation (7), G is the background signal in 
detected photons/sec/Hz per coherence area, M is the number of coherence 
areas, and D is the bandwidth in Hz of the optical filter. Maximization of 
the direct detection SNR with respect to the number of coherence areas occurs 
when M is selected so that the background and speckle terms are equal : 
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For this condition , the criterion for direct detection SNR to exceed the 
single-detector heterodyne large signal SNR of ^St is for >2T/GBD. 


SENSITIVITY ANALYSIS FOR DIAL MEASUREMENT USING 
HETERODYNE DETECTION WITH A SINGLE DETECTOR 

A comprehensive computer algorithm is used to calculate the expectation 
values of P^ . for various measurement conditions. Pressures, temperatures, 
and gas speciis densities are input from a midsummer midlatitude atmospheric 
model. Trace gas species densities can be modified using card inputs. Line 
absorption parameters are accessed from a comprehensive data base. Sources 
for line data are given in reference 6. At each altitude, molecular absorption 
at Vj and x >2 is calculated for each species by summing contributions from 
absorption lines in the vicinity of the laser frequency. Lorentz , Voigt, or 
Doppler line shapes are used at appropriate altitudes. Water vapor continuum 
absorption is added to the line absorption. Extinction due to particulate and 
molecular scattering is summed with molecular absorption at each altitude to 
give the total loss in each scattering cell. The integrated two-way loss is 
calculated by summing contributions from altitude layers between the laser and 
the cell being considered. The backscattering coefficient B is calculated 
by combining a Rayleigh term, which is small at infrared frequencies, and a Mb 
term. Mie backscatter and extinction for the cases shown here are calculated 
using parameters of Deirmendjian's Haze L size distribution (ref. 7). The 
vertical aerosol distribution is based on reference 8 with a ground level 
concentration (350 particles/cm^) corresponding to a 23-km visibility. 

Figure 1 shows the Mie contribution to the volume backscatter coefficient used 
in this study as a function of altitude for V = 927.61 cm - '*’. 

For all simulations presented, the overall system efficiency, which is 
the product of optical and quantum efficiency, is set at 12.5%. The system 
errors due to amplification and digitization are set at 0.1%. For aircraft 
and ground-based cases, telescope area is 0.1 m^ and bandwidth 10® Hz; for 
Shuttle cases, telescope area is 1 m^ and bandwidth 10^ Hz. Bandwidths have 
been constrained by atmospheric line parameters and by laser limits. Pulse 
energies are selected within practical laser constraints to match signal to 
bandwidth and the number of pulse pairs N is adjusted to obtain reasonable 
measurement accuracy. The "on" and "off" frequencies have been selected to be 
within the range of a single rare isotope multiatmospheric CO laser line. Us< 
of rare isotope lines minimizes interference by atmospheric c8^ and use of 
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closely spaced "on" and "off" frequencies minimizes error due to variations in 
backscatter and interfering species absorption. 


Figure 2 displays O 3 concentration in STP-ppb and measurement error in 
the same units versus altitude for measurement from Shuttle at 250 km. The 
pulse energy is 5 joules per pulse with 1000 pulse pairs per measurement. 
Telescope area is 1 m^ , bandwidth is 10^ Hz, resolution is 1.5 km below 30 km 
and 3 km above 30 km. Simulations are for an "off" frequency of 1058.01 cm - -*- 
and "on" frequencies of 1058.17, 1058.11, and 1058.19 cm - -*-. It should be 
noted that a fine tuning of the "on" frequency results in a sharp variation 
of altitude at which the best measurement can be accomplished. 

Figure 3 displays ozone concentration and error in STP-ppb for a measure- 
ment from aircraft at 10.5 km looking either upward or downward. Pulse energy 
is 0.05 joules, telescope area is 0.1 m 2 , bandwidth is 10® Hz, and "on" fre- 
quency is 1058.20 cm . In the case of aircraft measurements the "on" frequency 
selection is simplified since the signal does not have to pass through the ozone 
bulge. Two cases are presented in order to illustrate the strong influence of 
range resolution on number of pulses required for a measurement. The solid er- 
ror line is for 3000 pulse pairs with a resolution of 0.5 km below 6 km, 1 km 
from 6 to 10 km, and 3 km above 10 km. The dashed error line is for 100 pulse 
pairs with a resolution of 1.5 km. The errors are nearly identical at low al- 
titude, which is expected (see equation ( 6 )) since the number of pulses has been 
increased by approximately the cube of the inverse of the range resolution. 

Figure 4 displays NH 3 concentration and error in STP-ppb for 300 pulse 
pairs with operation from Shuttle at an altitude of 250 km with 1 joule per 
pulse and from aircraft at an altitude of 10 km with 0.05 joules per pulse. 

Since NH 3 is not heavily attenuated, only one set of frequencies is required. 

The resolution for both cases is 1.5 km. 


Figure 5 displays NH 3 concentration and error in STP-ppb for a ground- 
based measurement. The vertical NH 3 distribution is the same as shown in 
figure 4. For this case, the telescope area is 0.1 m 2 , bandwidth is 10® Hz, 
"on" frequency is 927.32 cm - -*-, and "off" frequency is 927.61 cm - -*-. The error 
curves are for a range resolution of 1.5 km using 300 pulse pairs. The two 
error curves are for laser energies of 0.05 J and 0.5 J. Note that for the 

/ 

cases in which E = 0.5 J, the large signal limit 


1: 


P i j + BT 


is achieved 


throughout the entire measurement range and therefore the error is nearly 
independent of altitude. 


Figure 6 is for the same conditions as shown in figure 5 with the excep- 
tion that the NH 3 concentration has been reduced by a factor of 5. The error 
is plotted only for the E = 0.05 J case. Comparison of figure 6 with the 
0.05 joule curve in figure 5 illustrates the effect of concentration on 
measurement error. At low altitude, P^j >BT and the measurement errors are 
independent of concentration. At higher altitudes P-j_j <BT and the increased 
integrated absorption at the higher concentration (fig. 5) results in an 
increased measurement error for that case. These NH 3 distributions are 
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typical values from Langley Research Center ground-based infrared heterodyne 
radiometer measurements (ref. 9) . 


Figure 7 displays percent error in measurement of water vapor concentra- 
tion versus altitude for measurement from 250 km for a midsummer, midlatitude 
water vapor distribution. Measurement conditions are 1 joule per pulse, 100 
pulse pairs, 10^ Hz bandwidth, and 1 m 2 telescope. Two "on" frequencies are 
used at 948.25 and 948.30 cm”-'- . The "off" frequency is at 948.35 cm - ^. Trop- 
ospheric water vapor is particularly amenable to measurement from space since 
the concentration rapidly increases towards the ground. Thus, absorption is 
high at low altitudes and the integrated absorption from space to low alti- 
tudes is relatively low. For measurement of water vapor at higher altitudes, 
the number of pulses would have to be increased. 


CONCLUDING REMARKS 

High vertical resolution measurement of atmospheric trace species can be 
achieved using CO^ laser DIAL with heterodyne detection. This study indicates 
that maximum sensitivity at minimum laser energy per measurement requires 
multiple pulse operation with the energy per pulse selected so that the measure' 
photon rate is approximately equal to the detector IF bandwidth. Measurement 
sensitivities can be maximized and interference effects minimized by fine ad- 
justment of measurement frequencies using the tunability of high-pressure laser 
The use of rare isotope lasers minimizes loss due to CO^ atmospheric absorption 
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Figure 1.- Aerosol volume backscattering coefficient as a function of 

altitude for V = 927.61 cm -1 . 



0 1 0 20 30 40 50 

ALTITUDE, km 

Ar = 1.5 km j Ar = 3 km 


Figure 2.- Ozone measurement from 250 km. 
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ALTITUDE IN KILOMETERS 

Figure 3.- Ozone measurement from 10.5 km. 



Ar = 1. 5 km 

Figure 4.- NH^ measurement from 10 and 250 km. 


-6 km 
-10 km 


565 




A r = 1.5 km 

Figure 5-- NH^ measurement from ground for E = 0.5 and 0.05 J. 



Figure 6.- NH^ measurement from ground for E = 0.05 J. 
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Figure 7.- Water vapor measurement from 250 km 





